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ABSTRACT 

Seitenzahl et al. have predicted that roughly three years after its explosion, the light we receive from a Type 
la supernova (SN la) will come mostly from reprocessing of electrons and X-rays emitted by the radioactive 
decay chain 57 Co -A 57 Fe, instead of positrons from the decay chain 56 Co -A 56 Fe that dominates the SN 
light at earlier times. Using the Hubble Space Telescope , we followed the light curve of the SN la SN 2012cg 
out to 1055 days after maximum light. Our measurements are consistent with the light curves predicted by the 
contribution of energy from the reprocessing of electrons and X-rays emitted by the decay of 57 Co, offering 
evidence that 57 Co is produced in SN la explosions. However, the data are also consistent with a light echo 
~ 14 mag fainter than SN 2012cg at peak. Assuming no light-echo contamination, the mass ratio of 57 Ni and 
56 Ni produced by the explosion, a strong constraint on any SN la explosion model, is 0.043 +q q}^, roughly 
twice Solar. In the context of current explosion models, this value favors a progenitor white dwarf with a mass 
near the Chandrasekhar limit. 

Subject headings: nuclear reactions, nucleosynthesis, abundances - supernovae: general - supernovae: indi¬ 
vidual (SN 2012cg) 


1. INTRODUCTION 

The declining light curves of Type la supernovae (SNe la), 
as well as some core-collapse SNe, are powered by the ra¬ 
dioactive decay of 56 Ni -A 56 Co -A 56 Fe, which releases 
y-rays that thermalize with the expanding ejecta, to be re¬ 
emitted in the optical and near-infrared (Truran et al. 1967; 
Colgate & McKee 1969). As the ejecta continue to expand, 
the column density decreases so that the ejecta become pro¬ 
gressively more transparent to y-rays. However, 19% of 56 Co 
decays proceed via positron emission, and after ~ 200 days 
these positrons become the dominant energy source and con¬ 
tinue to heat the ejecta (e.g., Arnett 1979; Chan & Lingenfel- 
ter 1993; Cappellaro et al. 1997; Milne et al. 1999). Assuming 
100% trapping of the positrons emitted by these decays, the 
observed light curves of SNe la 300-600 days after explo¬ 
sion (e.g., Sollerman et al. 2004; Lair et al. 2006) are con¬ 
sistent with the decay of 56 Co. However, the apparent slow¬ 
down of SN la light curves past 800 days, as seen in, e.g., 
SN 1992A (Cappellaro et al. 1997), SN 2003hv (Leloudas 
et al. 2009), and SN 201 lfe (Kerzendorf et al. 2014), indicates 
that even complete trapping of the 56 Co positrons cannot ex¬ 
plain the observed luminosity; additional radioactive heating 
channels may be required. Seitenzahl et al. (2009) suggested 
that the slower decays of 57 Co —>■ 57 Fe (7, /2 « 272 days) 
and 55 Fe —» 55 Mn (/j /2 ~ 1000 days) will produce internal- 
conversion electrons, Auger electrons, and X-rays that would 
largely be trapped, deposit their energy, and slow down the 
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decline of the light curve beginning ^ 800 days after maxi¬ 
mum light. 

While there is a wide consensus that SNe la are the result of 
a thermonuclear explosion of a carbon-oxygen white dwarf, 
it is still unclear how the white dwarf is detonated and con¬ 
sumed (for a recent review, see Maoz et al. 2014). Ropke et al. 
(2012) have used the leptonic energy injection mechanisms 
due to the decays of 57 Co and 55 Fe described above to provide 
a testable prediction for two popular SN la explosion scenar¬ 
ios: the “delayed-detonation” and “violent merger” models 
(Khokhlov 1991; Pakmor et al. 2011, 2012). In the first sce¬ 
nario, a white dwarf with a mass close to the Chandrasekhar 
limit (Mch ~ 1.4 M Q ) first goes through a subsonic defla¬ 
gration phase before transitioning to a supersonic detonation. 
In the second scenario, the white dwarf detonates due to a 
merger with a second, slightly lower-mass white dwarf, and 
both stars have sub-Chandra masses (1.1 and 0.9 M 0 ). Each 
of these scenarios predicts that different amounts of 57 Co and 
55 Fe will be created during the explosion. For example, in 
the merger scenario, where the densities of the white dwarfs 
are < 2 x 10 8 g cm 3 , most of the iron-group elements, in¬ 
cluding 55 Co (which then decays to the longer-lived 55 Fe), 
are produced in incomplete Si burning and a-rich freezeout 
from nuclear statistical equilibrium. In contrast, the higher 
density of the delayed-detonation white dwarf leads to “nor¬ 
mal” freezout (which sets in at ~ 2 x 10 8 g cm -3 ), which 
produces more 55 Co (Thielemann et al. 1986). The different 
abundances of 55 Fe and 57 Co then lead to bolometric light 
curves with appreciably different slopes at > 1200 days after 
maximum light. 

The late-time observations of the core-collapse SNe 
1998bw (Seitenzahl et al. 2009) and 1987A (Seitenzahl et al. 
2014) were shown to be consistent with the lep tonic heating 
mechanism suggested by Seitenzahl et al. (2009). When this 
paper was submitted, this prediction had not been tested for 
any normal SN la (e.g., SN 200lei, observed by Stritzinger 
& Sollerman 2007, as well as the SNe la collected by Milne 
et al. 2001 and Lair et al. 2006, only go out to ~ 600 days after 
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Fig. 1 HST color composite of NGC 4424 and SN 2012cg (in the center of the white box). The image is composed of all F160W (red), F814W (green), and 
F555W (blue) images taken as part of the GO-12880 program. The inset shows two blow-ups of the vicinity of SN 2012cg taken with the F350LP filter a year 
apart on 2013 December 24.9 (GO-12880) and 2014 December 16.5 (GO-13799). Each panel is 2" on a side. The chip gap of the WFC3/UVIS detector has 
been filled in based on data from the F160W images. North is up and East is to the left. 


explosion). Tantalizingly, the very last data points on the light 
curves of SN 1992A 950 days; Cappellaro et al. 1997) and 

SN 2003hv 700 days; Leloudas et al. 2009) might be show¬ 

ing the beginning of a slow-down of their optical light curves. 
Very recently, Fransson & Jerkstrand (2015) published their 
analysis of the spectrum of SN 201 lfe taken 1034 days after 
explosion by Taubenberger et al. (2015). At this epoch, they 
find strong evidence of the need for energy injection by 57 Co 
and, in order to reproduce the observed flux level of the spec¬ 
trum, require a production ratio of 57 Ni/ 56 Ni that is roughly 
2.8 times Solar. Here, we test the predictions made by Seiten- 
zahl et al. (2009) and Ropke et al. (2012) by following the 
late-time lightcurve of the SN la SN 2012cg. 

SN 2012cg was discovered on 2012 May 17 (UT) in the 
nearby spiral galaxy NGC 4424 (at a stellar-kinematics based 
distance of 15.2 ± 1.9 Mpc; Cortes et al. 2008) by the Lick 
Observatory Supernova Search (Filippenko et al. 2001; Kan- 
drashoff et al. 2012) and spectroscopically classified as a SN 
la by Cenko et al. (2012) and Marion et al. (2012). Hubble 
Space Telescope (HST) pre-explosion images of NGC 4424 
revealed no source at the location of the SN down to limits of 
My ~ —6.0 and Mj ~ —5.4 mags and excluded most super¬ 
giants as potential binary companions (Graur & Maoz 2012). 
Chomiuk et al. (2012) did not detect the SN in the radio and 
constrained any possible circumstellar material to either lie at 
distances > 10 16 cm or be distributed in thin shells < 10 15 cm 
wide. 


2. OBSERVATIONS 

We imaged SN 2012cg (a = 12 h 27 m 12.83 s , 8 = 

+09°25'13.1") with the HST Wide Field Camera 3 (WFC3) 
wide-band filters F350LP, F555W, F814W, and F160W un- 
der HST program GO-12880 (PI Riess) on 2013 December 
24.9; 2014 January 1.2, 8.8, 17.1, 23.6, and 28.9; February 
4.0, 11.8, 20.8, and 27.5; and March 9.6. At these dates the 
SN was ~ 570-640 days past maximum light. We subse¬ 
quently observed SN 2012cg under HST program GO-13799 
(PI Graur) with the WFC3 filter F350LP on 2014 December 
16.5, 2015 February 6.9, and 2015 April 26.6, when the SN 
was « 925-1055 days beyond peak light. Here, we concen¬ 
trate on observations carried out with the long-pass F350LP 
filter. Centered at wavelength Ao = 584.6 nm and with a width 
of 475.8 nm, F350LP is the second-widest WFC3 filter, after 
F200LP , covering a wavelength range of ~ 300-1000 nm. 

In Figure 1, we show a color composite of NGC 4424 and 
SN 2012cg along with an inset that shows F350LP observa¬ 
tions of the SN taken a year apart. While the SN has faded 
appreciably between the two visits, it is still visible > 900 
days after explosion. We present our measurements of the 
light curve of SN 2012cg in Figure 2. 

3. ANALYSIS 

We performed aperture photometry of SN 2012cg in a 0.4"- 
diameter aperture using the hstsnphot 7 program used by the 

7 https ://github.com/srodney/hstsntools 





Discovery of 57 Co in SN 2012cg 


3 


TABLE 1 

Late-time photometry of SN 2012cg 


Date 

MJD 

(days) 

Phase 

(days) 

Filter 

Exposure Time 
(s) 

Magnitude 
(AB mag) 

A(Magnitude) 
(AB mag) 

Flux 

(MJy) 

A(Flux) 

(MJy) 

HST Program 

2013 Dec 24.9 

56650.9 

567.9 

F350LP 

1200 

23.094 

0.039 

2.10 

0.08 

12880 

2013 Dec 24.9 

56650.9 

567.9 

F555W 

936 

22.750 

0.019 

2.88 

0.05 

12880 

2014 Jan 01.2 

56658.2 

575.2 

F350LP 

1200 

23.217 

0.040 

1.88 

0.07 

12880 

2014 Jan 01.2 

56658.2 

575.2 

F814W 

900 

22.867 

0.034 

2.59 

0.08 

12880 

2014 Jan 08.8 

56665.8 

582.8 

F350LP 

1260 

23.271 

0.041 

1.78 

0.07 

12880 

2014 Jan 08.8 

56665.8 

582.8 

F160W 

906 

21.35 

0.14 

10.5 

1.3 

12880 

2014 Jan 17.1 

56674.1 

591.1 

F350LP 

1200 

23.341 

0.042 

1.67 

0.06 

12880 

2014 Jan 17.1 

56674.1 

591.1 

F160W 

906 

21.51 

0.17 

9.0 

1.4 

12880 

2014 Jan 23.6 

56680.6 

597.6 

F350LP 

1230 

23.404 

0.043 

1.58 

0.06 

12880 

2014 Jan 23.6 

56680.6 

597.6 

F814W 

900 

23.173 

0.046 

1.95 

0.08 

12880 

2014 Jan 28.9 

56685.9 

602.9 

F350LP 

1200 

23.439 

0.043 

1.53 

0.06 

12880 

2014 Jan 28.9 

56685.9 

602.9 

F814W 

938 

23.291 

0.050 

1.75 

0.08 

12880 

2014 Feb 04.0 

56692.0 

609.0 

F350LP 

1260 

23.506 

0.044 

1.44 

0.06 

12880 

2014 Feb 04.0 

56692.0 

609.0 

F160W 

906 

21.80 

0.21 

6.9 

1.3 

12880 

2014 Feb 11.8 

56699.8 

616.8 

F350LP 

1260 

23.570 

0.045 

1.35 

0.06 

12880 

2014 Feb 11.8 

56699.8 

616.8 

F160W 

906 

21.79 

0.21 

7.0 

1.4 

12880 

2014 Feb 20.8 

56708.8 

625.8 

F350LP 

1230 

23.640 

0.047 

1.27 

0.06 

12880 

2014 Feb 20.8 

56708.8 

625.8 

F555W 

900 

23.132 

0.026 

2.03 

0.05 

12880 

2014 Feb 27.5 

56715.5 

632.5 

F350LP 

1230 

23.738 

0.048 

1.16 

0.05 

12880 

2014 Feb 27.5 

56715.5 

632.5 

F814W 

900 

23.566 

0.063 

1.36 

0.08 

12880 

2014 Mar 09.6 

56725.6 

642.6 

F350LP 

1260 

23.785 

0.050 

1.11 

0.05 

12880 

2014 Mar 09.6 

56725.6 

642.6 

F160W 

906 

21.96 

0.24 

5.9 

1.3 

12880 

2014 Dec 16.5 

57007.5 

924.5 

F350LP 

2688 

25.76 

0.22 

0.18 

0.04 

13799 

2015 Feb 06.9 

57059.9 

976.9 

F350LP 

2688 

25.94 

0.26 

0.15 

0.04 

13799 

2015 Apr 26.6 

57138.6 

1055.6 

F350LP 

2688 

26.29 

0.29 

0.11 

0.03 

13799 


Note. All photometry is measured in an aperture with a diameter of 0.4" (except for the final measurement on 2015 26.6, for which we 
used an aperture with a diameter of 0.24"). Phases are relative to 5-band maximum light on 2012 June 4.5. The measurements here have 
not been corrected for reddening by either Galactic or host-galaxy dust. 


CLASH and CANDELS SN surveys (Graur et al. 2014a; Rod¬ 
ney et al. 2014), which is based on the IRAF 8 apphot routine 
(Tody 1986). Due to the faintness of the SN on the last ob¬ 
servation epoch, 2015 April 26.6, we used a smaller aperture 
with a diameter of 0.24", along with a larger aperture correc¬ 
tion, to capture the low flux counts. Zero points and aperture 
corrections were measured by the Space Telescope Science 
Institute. 9 Our photometry of SN 2012cg, uncorrected for 
Galactic or host-galaxy reddening, is reported in Table 1 . 

We correct our measurements for both Galactic extinction 
along the line of sight to SN 2012cg and the extinction caused 
by dust in the SN host galaxy. The Galactic line-of-sight 
extinction toward SN 2012cg in F350LP , F555W , F814W , 
and F160W is 0.054, 0.059, 0.032, and 0.011 mag, respec¬ 
tively (Schlafly & Finkbeiner 2011). By using the Cardelli 
et al. (1989) reddening law, an E {B — V) = 0.18 mag redden¬ 
ing (Silverman et al. 2012; Marion et al. 2015), and assum¬ 
ing Ry = 3.1, we calculate host-galaxy extinctions of 0.515, 
0.582, 0.333, and 0.114 mag in the same filters. 

SN 2012cg peaked in the B band sometime during 2012 
June 2-5, at which point it reached an absolute magnitude 
similar to that of SN 201 lfe (as measured by Munari et al. 
2013). While it is similar to SN 201 lfe in many respects 
(e.g., in its color evolution around peak, as well as the ap- 
perance of its nebular spectra; see Amanullah et al. 2015, 
but cf. Marion et al. 2015 for a possible early excess of 
blue light. For more information on SN 201 lfe, see, e.g., 
Li et al. 2011; Chomiuk 2013; Kasen & Nugent 2013; Graur 

8 http://iraf.noao.edu/ 

9 http://www.stsci.edu/hst/wfc3/phot_zp_lbn 


et al. 2014b), it has a slightly broader light curve. Specifi¬ 
cally, by fitting early-time data, Silverman et al. (2012) mea¬ 
sured q ax = 2012 June 2.0 =b 0.75, M B = -19.73 ± 0.30 mag, 
and a ZLband light curve that was slightly narrower than SN 
201 lfe (the only work to reach such a conclusion). Munari 
et al. (2013) measured t^ ax = 2012 June 4.5, M B = —19.55 
mag, Am\ 5 :B = 1.039 (compared to 1.108 for SN 201 lfe), and 
a reconstructed bolometric light curve that reaches a similar 
(though slightly higher) peak luminosity as SN 201 lfe, with 
a slightly broader shape overall (see their figure 4). Marion 
et al. (2015) measured an earlier t% ax = 2012 June 3.3 =b 0.5, 
M b = —19.62 ±0.08 mag, and A m\$ )B = 0.86 ±0.02. Finally, 
Amanullah et al. (2015) found that t^ ax = 2012 June 2.6 ± 0.3 
and s'# = 1.12 ± 0.02, where the latter is the Perlmutter et al. 
(1997) stretch parameter. Here, we adopt the value of peak B- 
band light derived by Munari et al. (2013), so that all phases 
are reported relative to 2012 June 4.5. 

The decay of 56 Co to 56 Fe has a half-life time of = 
77.2 days (Junde et al. 2011). Assuming that the shape of 
the light curve is set by 56 Co decay, we fit an exponential 
function to the F350LP measurements at 570 < / < 640 days 
and find a 56 Co half-life of 82.7 ±5.3 days with a x 2 value of 
1.7/9 for nine degrees of freedom. A similar fit to the F160W 
data results in a half-life of 69 ±26 days with x 2 = 0.6/3. 
These values are consistent with the half-life time obtained 
in the laboratory, and with each other. However, they also 
hint that at > 500 days, we may be seeing a redistribution 
of light from the infrared to the optical, as previously seen 
in observations of SN 1992A and SN 2003hv (though based 
on only one and two measurements, respectively; Cappellaro 
et al. 1997; Leloudas et al. 2009). 
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In Figure 2, we extrapolate the F350LP fit out to 1100 days 
to represent the expected light curve due solely to the decay 
of 56 Co to 56 Fe. The three F350LP measurements at t > 900 
days are brighter than the extrapolated fit by factors of 1.8- 
3.2. To ascertain at what significance these measurements are 
inconsistent with the null hypothesis that the light curve is due 
solely to 56 Co decay, we compare between the three t > 900 
days measurements and the model, while taking into account 
the uncertainties of both. When fitting these measurements, 
the values for the index and scaling of the decaying exponen¬ 
tial function that resulted from fitting the 570 < t < 640 days 
measurements receive a x 2 value of 9.9/1, which calls for a 
rejection of the null hypothesis at a significance of > 99.8% 
(i.e., a p -value of ~ 0.002, comparable to a 3d detection). 

Next, we test whether the t > 900 days F350LP measure¬ 
ments are consistent with the predictions for the delayed- 
detonation and violent merger progenitor scenarios. The orig¬ 
inal predictions are bolometric light curves of SN 201 lfe, 
while we measure the optical component of the light curve 
of SN 2012cg captured by the F350LP filter. Although the 
predictions were made for a different SN la, SN 2012cg was 
observed to have a similar intrinsic maximal brightness as 
SN 201 lfe but a slightly broader light curve. We solve both 
of these potential difficulties by comparing between the pre¬ 
dicted bolometric light curve due to 56 Co decay for SN 201 lfe 
and our fit to the 570 < t < 640 days F350LP observations of 
SN 2012cg and scaling down the predicted bolometric light 
curves accordingly. The F350LP measurements are consistent 
with the resultant scaled-down predictions but are not precise 
enough to discriminate with certainty between the two pro¬ 
genitor scenarios. 

Here, we have assumed that the emission in F350LP scales 
with the bolometric luminosity. Milne et al. (2001) have 
shown that at 50 < t < 600 days, a constant fraction of the 
luminosity is emitted in the V band, which dominates the 
3500-9700 A wavelength range (the same range spanned by 
F350LP). Past ~ 500 days, the V band becomes dominated 
by Fe lines (see Taubenberger et al. 2015), which, in the ab¬ 
sence of the infrared catastrophe (cf. Fransson & Jerkstrand 
2015), are the primary coolant of the SN ejecta and hence also 
proportional to the total luminosity of the SN. 

As noted above, the two progenitor scenarios tested here 
differ, among other things, in the amounts of radioactive iron- 
group elements produced during the explosion. Thus, a mea¬ 
surement of the ratio between the amount of 56 Ni and 57 Ni 
produced during the explosion would provide a strong con¬ 
straint for any progenitor and explosion model. Here, we 
compute the luminosity contributed by the decays of 56 Co, 
57 Co, and 55 Fe by using the solution to the Bateman equation 
in the following form (following Seitenzahl et al. 2014): 

= 2 - 221 ?d^i^sv exp( -^' ) x 1043 erg s_ ‘- 

( 1 ) 

where B = 0.235 is the scaling factor described above; A is the 
atomic number of the decaying nucleus; Xa is the inverse of 
the half-life time of the decay chain; qA is the average energy 
per decay carried by charged leptons and X-rays; t is the time 
since explosion; and we fit for the masses, M(A), of 56 Co 
and 57 Co. Our late-time measurements are not precise enough 
to also fit for the mass of 55 Fe, so we set its mass by using 
a ratio of M( 57 Co)/M( 55 Fe) « 0.8 (model rpc32; Ohlmann 



Fig. 2.— Light curve of SN 2012cg at 570 < t < 1055 days after 5-band 
maximum light. F350LP, F555W, F814W, and F160W luminosities appear 
as white, green, red, and magenta squares, as noted. The solid black curve 
is the best-fitting exponential function to the F350LP measurements in the 
range 570 < t < 640 days. The shaded region around this curve represents 
the confidence region of the fit given the 68% statistical uncertainties of the 
fitted parameters. The dashed purple and dotted-dashed orange curves are the 
predicted bolometric luminosities of SN 201 lfe for the delayed-detonation 
and violent merger progenitor scenarios, respectively, scaled down so that the 
predicted 56 Co-only decay matches the fit to the 570 < t < 640 days F350LP 
data. The F350LP measurements at t > 900 days indicate that the light curve 
is no longer dominated solely by the decay of 56 Co; instead, its decline is 
being slowed down by light emitted from the decay of 57 Co to 57 Fe. 

et al. 2014). The values of Xa and qA used here are enumerated 
in table 1 of Seitenzahl et al. (2009) and table 2 of Seitenzahl 
et al. (2014). For 57 Co, we also take into account the energy 
emitted by the decay through the 14.4 keV line (1.32 keV per 
decay). 

Several assumptions go into Equation 1 : 

1. y-rays are essentially free-streaming and their contribu¬ 
tion to the heating is negligible. 

2. Positrons, electrons, and X-rays are fully trapped and 
instantaneously deposit all their energy. 

3. The luminosity is equal to the rate of energy deposi¬ 
tion at every time, i.e., we do not take “freeze-out” into 
account (see, e.g., Kozma & Fransson 1998). 

While these assumptions clearly constitute approximations, 
we presume they are applicable and justified since: 

1. The dominant opacity to y-rays is Compton scattering. 
For homologous expansion, the column density, and 
therefore the opacity to y-rays, decreases as t~ 2 . Since, 
e.g., Cappellaro et al. (1997) find that for a generic SN 
la model the trapping fraction of y-rays at t = 200 days 
is already down to 0.5%, we conclude that the contribu¬ 
tion of y-rays is probably safely negligible past ~ 600 
days (with the exception of the very low-energy 14.4 
keV y-ray produced by 57 Co decay, which we include 
in our analysis). 

2. The calculations of Seitenzahl et al. (2015) demonstrate 
that at ~ 1000 days full trapping of X-rays produced 
in radioactive decay of Fe-group nuclei should be a 
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good approximation. The possible escape of positrons 
is more uncertain. However, at present there is no real 
indication for positron escape in normal SNe la. On 
the contrary, recent analyses of very late SN la obser¬ 
vations support full trapping of positrons even past 800 
days (see, e.g., Leloudas et al. 2009; Kerzendorf et al. 
2014; Fransson & Jerkstrand 2015). If positrons are 
fully trapped, the much lower-energy electrons would 
be trapped as well. 

3. Fransson & Jerkstrand (2015) demonstrated that, for 
the case of the similar SN 201 lfe, freeze-out effects 
are essentially negligible out to ~ 900 days in the V 
band, which dominates the emission in F350LP , tak¬ 
ing the Taubenberger et al. (2015) SN 201 lfe spectrum 
at 1034 days past explosion as our best guess for the 
spectral-energy distribution. Furthermore, their freeze- 
out corrections remain smaller than our measurement 
uncertainties even at ~ 1050 days. Freeze-out effects, 
therefore, constitute a small correction that would not 
change our conclusions qualitatively. 

With a x 1 value of 2.1 for 12 degrees of freedom, we find 
a best-fitting 56 Co mass of M( 56 Co) « 0.7 M 0 and a mass ra¬ 
tio of M( 57 Co)/M( 56 Co) = 0.043!qq} 7 . We show the lumi¬ 
nosity contributions from each decay chain in Figure 3. Re¬ 
stricting the fit to account only for 56 Co and 57 Co results in 
the same 56 Co mass but a mass ratio of M( 57 Co) /M( 56 Co) = 
0.068+QQ}g, larger by a factor of ~ 1.6 (with a % 2 value of 

2.5/12). We do not quote an uncertainty for the 56 Co mass 
as its value is degenerate with the scaling factor B. Without 
a full bolometric light curve, we cannot measure this quantity 
precisely, but note that it is consistent with 56 Ni masses mea¬ 
sured in other SNe la (Scalzo et al. 2014). The mass ratio, 
however, is independent of the scaling factor. In Section 4.1, 
below, we examine the value of the mass ratio in the context 
of expected yields from various SN la explosion models. 

Here, we have used Ry = 3.1 to correct for the host-galaxy 
reddening suffered by SN 2012cg. However, based on the ob¬ 
served colors of SN 2012cg around maximum light, Amanul- 
lah et al. (2015) find a lower value of Ry = 2 . 7 +q 7. While this 
value is consistent with the one we use here, we note that the 
value of Ry has no effect on the value of M( 57 Co) /M( 56 Co), 
as it depends only on the slope of the F350LP photometry, not 
their absolute values. 

4. DISCUSSION 

In this section, we first place the M( 57 Co)/M( 56 Co) mass 
ratio found above in the context of current SN la explosion 
models. Then, we address two other possible explanations for 
the slow-down of the decline of the light curve at > 900 days: 
contamination by light echoes or the light from a surviving 
binary companion. We show the effects of these possible con¬ 
taminants in Figure 4. 

4.1. Expected 57 Co yields from SN la explosion models 

Here, we assume that nearly all of the mass of 
the iron-group elements produced by the explosion orig¬ 
inates in the decay of various Ni isotopes, so we 
can approximate M( 57 Ni)/M( 56 Ni) = M( 57 Co)/M( 56 Co) = 
M( 57 Fe)/M( 56 Fe) (Truran et al. 1967; Woosley et al. 1973). 
Under this approximation, the mass ratio we measure here is 



Fig. 3.— Luminosity contributions from the decays of 56 Co (blue dashed), 
57 Co (red dotted-dashed), and 55 Fe (green dotted). The total luminosity pro¬ 
duced by these decay chains (black solid) fits the F350LP measurements with 
a x 2 value of 2.1 for 12 degrees of freedom. 

roughly twice the corresponding 57 Fe/ 56 Fe ratio measured for 
the Sun (Asplund et al. 2009). This value is slightly higher 
than the predictions of near-Mch explosion models: W7 of 
Iwamoto et al. (1999) predicts 1.7 times Solar, the pure tur¬ 
bulent deflagration models of Fink et al. (2014) predict ~ 1.5 
times Solar, the delayed-detonation models of CO WDs of 
Seitenzahl et al. (2013) predict ~ 1.3 times Solar, the delayed- 
detonation models of WDs with carbon-depleted cores of 
Ohlmann et al. (2014) predict ~ 1.1 times Solar, and the grav¬ 
itationally confined detonation models of Meakin et al. (2009) 
predict 0.8-1.2 times Solar. 

There is a trend that models with larger progenitor metal- 
licity (see Seitenzahl et al. 2013) and models that experi¬ 
ence enhanced in-situ neutronization (e.g., W7) yield higher 
57 Ni/ 56 Ni ratios. This is naturally explained by the fact that 
57 Ni is a neutron-rich isotope (28 protons, 29 neutrons) and 
its production relative to the self-conjugate 56 Ni (28 protons, 
28 neutrons) therefore tends to be enhanced in neutron-rich 
environments. Model yields where the primary WD is sig¬ 
nificantly less massive than Mq h predict even lower 57 Ni/ 56 Ni 
ratios: the merger model of a 1.1 M 0 CO WD with a 0.9 M Q 
CO WD of Pakmor et al. (2012) predicts ~ 1.1 times Solar 
and the pure detonation models of CO and ONe WDs of Mar- 
quardt et al. (2015) predict around 0.3-0.5 times Solar. 

Parikh et al. (2013) have found that the 57 Ni/ 56 Ni pro¬ 
duction ratio is very robust (at the 10% level) to individual 
changes of nuclear reaction rates by a factor of 10; rate un¬ 
certainties are therefore unlikely to solely account for the dif¬ 
ference. Therefore, the high 57 Co/ 56 Co ratio we find appears 
to point toward a near-Mch progenitor. This is in agreement 
with Marion et al. (2015), who interpret the observed early 
UV excess of SN 2012cg as evidence that the SN occurred 
as an explosion of a near-Mch WD in the single-degenerate 
channel. 

4.2. Contamination by a light echo? 

Light echoes have been found for other SNe la at similar 
phases hundreds of days after their explosion, e.g., SN1991T 
(Schmidt et al. 1994; Sparks et al. 1999), SN1995E (Quinn 
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et al. 2006), SN1998bu (Cappellaro et al. 2001), SN2006X 
(Wang et al. 2008), SN2007af (Drozdov et al. 2015), and 
SN2014J (Crotts 2015). Typically, these light echoes were be¬ 
tween 10-12 mag fainter than the SNe at peak. Due to obser¬ 
vational bias, these observed light echoes are at the bright tail 
of the distribution of expected light-echo brightnesses (Patat 
2005). Therefore, the excess flux of SN 2012cg at > 900 days, 
which is ~ 14 mag fainter than the SN at peak, is consistent 
with a light echo. In the following, we investigate this possi¬ 
bility in more detail. 

There is no evidence of a resolved light echo at the loca¬ 
tion of SN 2012cg in the three > 900-day epochs. Using 
Dolphot , 10 a modified version of HSTphot (Dolphin 2000), 
we subtracted a model of the F350LP point-spread function 
(PSF) from the location of the SN. In each epoch, a visual 
inspection of the residual map produced by Dolphot revealed 
no resolved features. To quantify this statement, we compared 
the root-mean-square (RMS) of the pixel values at the location 
of the PSF-subtracted SN, measured in rectangular apertures 
of 7 x 7, 5 x 5, and 3x3 pixels 2 , to the distribution of RMS 
values measured in 100 apertures at random locations up to 
15 pixels away in any direction from the SN. In each of the 
three epochs, the RMS value was within the 16th and 84th 
percentiles of the distribution, i.e., consistent with the back¬ 
ground noise. 

We now address the possibility of a non-resolved light echo 
within the PSF. Let us first consider the shape of the light 
curve. As we have argued above, the observed excess flux 
surface brightness is well within the range expected for light 
echoes. With the reasonable assumption that the scattering 
dust structures do not vary significantly on the length scales 
associated with the light echo, the surface brightness of light 
echoes decays as t~ 2 (e.g., Chevalier 1986; Sugerman 2003). 
As the light echo is not resolved, we must also take into ac¬ 
count that its angular size, which is proportional to the cir¬ 
cumference, increases with time, so that the light echo should 
fade as t~ l . We test this hypothesis by fitting a combination of 
56 Co decay and a light echo decaying as t~ l to the data, with 
the light echo normalization set as a free parameter. Figure 4 
shows that this model, with x 2 = 1 . 8 / 12 , fits the data well. 

Next, we consider whether the observed colors are consis¬ 
tent with the light echo above. Light-echo flux is dominated 
by the light of the SN around peak, so to first order the color 
of the light echo is similar to the color of the SN at peak. 
In addition, the scattering by the dust favors blue light, which 
can shift the color of the light curve bluewards by a few tenths 
of a magnitude, depending on the scattering angle (e.g., Rest 
et al. 2012 ). 

At 0.3 days, Amanullah et al. (2015) measure V — i = 
—0.220 ±0.045 mag (Vega, corrected for Galactic extinction). 
Our photometry allows us to measure the F555W—F814W 
color (similar to V — I) of SN 2012cg at 570-640 days. As 
our F555W and F814W measurements were not taken on the 
same days, we use the python AstroML package (Ivezic et al. 
2014) 11 to perform Gaussian process regression and measure 
the color at any phase within this range. At 600 days, we mea¬ 
sure F555W — F814W = 0.140 d= 0.055 mag (Vega, corrected 
for Galactic extinction). Hence, at 600 days, SN 2012cg was 
redder than it was at peak. However, at 600 days there is 
no excess flux, assuming the light curve is powered solely 

10 http://americano.dolphinsim.com/dolphot/ 

11 http://www.astroml.org/ 



Fig. 4.— Possible contaminants of the light curve of SN 2012cg. The ex¬ 
pected luminosity range of a bright light echo is shown as a green patch, while 
56 Co decay (blue dashed) and a faint light echo that decays as t~ l (orange 
dotted-dashed) can combine (black solid) to fit the F350LP measurements 
(white squares). The dotted red curve shows the maximal contribution to the 
SN light from a red-giant or ~ 7 M© main-sequence binary companion, ac¬ 
cording to 2 a upper limits on the luminosity of such a companion measured 
from HST pre-explosion images. 

by 56 Co decay (see Figure 2). The combination of 56 Co de¬ 
cay and a faint light echo shown in Figure 4 would mean a 
flux excess of 11%, relative to the 56 Co component. There¬ 
fore, any contamination by a faint light echo would only min¬ 
imally change the observed combined color. Thus, the ob¬ 
served color of SN 2012cg at 600 days is consistent with ei¬ 
ther no light echo or a light echo faint enough that it leaves no 
trace at this phase. At > 900 days, where the excess flux dom¬ 
inates the observed flux and the light-echo hypothesis could 
be tested best with the color, we only have observations in 
one filter. 

Without much fine-tuning, a light echo can explain the ob¬ 
served colors and light-curve shape. Thus, with no spectrum 
or color information at late times, we cannot conclusively rule 
out contamination of our measurements by a light echo. It is 
a coincidence of nature that the functional forms of the light- 
echo and 57 Co-decay luminosities exhibit the same behavior 
in the time range probed here, as can be seen by comparing 
Figures 3 and 4. Unfortunately, further observations of SN 
2012cg cannot inform this issue. Our final F350LP measure¬ 
ment already had a signal-to-noise ratio of 3.3; by now, the 
SN should be as faint or fainter than the background light of 
NGC 4424. 

4.3. Contamination by a surviving binary companion? 

To test whether our measurements at / > 900 days indicate 
that the SN had faded to the point where a surviving binary 
companion could shine through, we first fit all F350LP mea¬ 
surements with a combination of an exponential function, for 
the SN light curve, and a constant luminosity from the pre¬ 
sumptive companion. The best-fitting companion has a lumi¬ 
nosity of (4.4^j 4 ) x 10 3 Lq ( x 2 = 2.4/12), consistent with a 
red giant. However, pre-explosion HST Wide-Field Planetary 
Camera 2 images of the location of the SN have produced 
2cr upper limits of F606W = 25.5 and F814W = 25.8 mag 
(Graur & Maoz 2012). Taking into account Galactic and host- 
galaxy extinctions for these filters results in a 2(7 upper limit 
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on the F350LP luminosity of a companion star of 840 L 0 , 
as shown in Figure 4. Such a luminosity is consistent with 
either a red giant or a ~ 7 M 0 main-sequence star, but con¬ 
stitutes only ~ 10% of our measured F350LP luminosities at 
t > 900 days. Our upper limit is consistent with the conclu¬ 
sion of Marion et al. (2015) that the pre-maximum UV excess 
observed in SN 2012cg is the possible signature of interaction 
between the SN ejecta and a ^ 6 M 0 main-sequence compan¬ 
ion. We conclude that, while we cannot rule out the presence 
of a surviving companion star, even if it is present, its contri¬ 
bution to the light curve is still negligible. 


5. CONCLUSIONS 

In conclusion, observations of the SN la SN 2012cg more 
than 900 days after its explosion allow us to reject, at a sig¬ 
nificance of >99.8% (i.e., 3a), the null hypothesis that the 
light curve of this SN la is powered solely by the radioactive 
decay chain 56 Ni -A 56 Co -A 56 Fe. Instead, we have shown 
that our measurements are consistent with predictions that at 
these late times, the light curve becomes dominated by repro¬ 
cessing of internal-conversion and Auger electrons, as well as 
X-rays, emitted by the decay of 57 Co into 57 Fe, as predicted 
by Seitenzahl et al. (2009). This constitutes the first evidence 
for the synthesis of radioactive 57 Co during a SN la explosion. 

We have tested whether the late-time behavior of the light 
curve of SN 2012cg could be due to contamination by either 
a light echo or the luminosity from a surviving binary com¬ 
panion. Pre-explosion imaging sets an upper limit on the lat¬ 
ter, which makes it a negligible source of contamination. A 
combination of 56 Co decay and a faint light echo, ~ 14 mag 
fainter than SN 2012cg at peak, fits the data well. Unfortu¬ 
nately, without spectra or colors at late times, we cannot rule 
out the possibility of light-echo contamination. 

Assuming no contamination by light echoes, we have mea¬ 
sured a ratio of 0.043 +qq} 7 between the masses of 57 Co and 
56 Co produced by the decay of Ni isotopes synthesized in the 
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